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Abstract 

The B*Btt coupling is a fundamental parameter of chiral effective Lagrangian 
with heavy-light mesons and can constrain the chiral behavior of Bb and the 
B — > nlv form factor in the soft pion limit. We compute the B*Bn coupling with 
the static heavy quark and the 0(a)-improved Wilson light quark. Simulations 
are carried out with n/ = 2 unquenched 12 3 x 24 lattices at j3 = 1.80 and 16 3 x 32 
lattices at (5 = 1.95 generated by CP-PACS collaboration. To improve the statistical 
accuracy, we employ the all-to-all propagator technique and the static quark action 
with smeared temporal link variables following the quenched study by Negishi et 
al. These methods successfully work also on unquenched lattices, and determine 
the B*Btt coupling with 1-2% statistical accuracy on each lattice spacing. 
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1 Introduction 



One of the major subjects in particle physics is to determine the CKM matrix elements 
in order to test the standard model and find a clue to the physics beyond. While the 
precision of the experimental data from B factories having been improving significantly, 
there are still large uncertainties in the CKM matrix elements due to the theoretical 
errors, which includes those in the lattice determination of the weak matrix elements for 
the B mesons. 

It is often the case that a symmetry helps to obtain nonperturbative results in field 
theories. For example, the chiral Lagrangian based on the approximate chiral symmetry 
can help to understand the quark mass dependence of the light mesons and also to derive 
nontrivial relations between different physical quantities related by the chiral symmetry. 
For the B mesons, there is another symmetry called 'heavy quark symmetry' which ap- 
pears in the limit of infinitely large quark mass. Based on this symmetry one can construct 
the heavy meson effective theory, which gives a systematic description of the heavy-light 
mesons including 1/M corrections. Using this effective theory one can understand the 
light quark mass dependence of various physical observables of the B meson weak matrix 
elements and can also derive nontrivial relations between different quantities, provided 
the low energy constants being determined from some method. 

The heavy meson effective Lagrangian has single low energy constant at the leading 
order of the 1/M expansion. This constant, is called the B*Btc coupling. Once 
the B*Btt coupling is determined, the heavy meson effective theory can predict various 
quantities which are important for CKM phenomenology [1] . For example the light quark 
mass dependence of the B meson decay constant and the bag parameter can be determined 
as 

/ 3 m 2 \ 

f Bd = F\l + -(l + 3g 2 b )^-^\og(m 2 JA 2 )j + analytic terms, (1) 

f 3 m 2 \ 

B Bd = B (l + -(1 - logK/A 2 ) J + analytic terms. (2) 

F and B are the low energy constants associated with these operators, and correspond to 
those quantities in the chiral limit of the light quark. Also the form factor f + (q 2 ) for the 
semileptonic decay B — > nlu can be expressed in terms of the B* meson decay constant 
fs* and cjb as 



f + {q 2 



f. 



B* 



2/. 



m B * _ m B * \ f B 
v-k-A m B J f B * 



(3) 



where v is the velocity of the B meson, k is the pion momentum, and A = m B * — m B . 
Therefore it is quite important to determine the B*Bir coupling very precisely from lattice 
QCD simulations. For this purpose, one of the promising approaches is to use the heavy 
quark effective theory (HQET) with nonperturbative accuracy including 1/M corrections. 
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HQET allows systematic treatment of the b quark in the continuum theory where 1/M 
corrections can also be systematically included with nonperturbative accuracy. 

Despite its usefulness, it is very difficult in practice to calculate the matrix elements 
for heavy-light systems with HQET [21 El E] • This is because in the heavy-light system 
the self-energy correction to the static quark gives a significant contribution to the energy, 
which results in an exponential growth in time of the noise to signal ratio of the heavy-light 
meson correlators. In fact, recent results of g^ are 

3oo = 0.51±0.03 stat ±0.11 sys for n f = % (4) 

= 0.51 ± 0.10 s tat for n f = 2 % (5) 

which have about 5% and 15% statistical errors for quenched and unquenched cases, re- 
spectively. An alternative method which extracts g^ from the B quark potential was 
proposed in Ref. [5], but such accuracies would not be sufficient to test new physics. 
Therefore significant improvements for statistical precision in HQET are needed. Fortu- 
nately the two techniques to reduce the statistical error are developed recently, which are 
the new HQET action [HI [7] with HYP smearing [8] and the all-to-all propagators [9] with 
the low mode averaging [TQl E] ■ Negishi et al. p2] tested applicability of these methods 
on a quenched lattice, and found that the statistical accuracy is drastically improved as 

<?oo = 0.517(16) sta t. for n/ = 0, (6) 

namely at 2% level, even with a modest number of configurations. 

Our final goal is to extend the above strategy to unquenched simulations and give 
a precise value of the B*Bn coupling g^ with 2 + 1 flavors in the continuum limit. In 
this paper, we study the static B*Bn coupling in nf = 2 unquenched QCD combining 
two techniques of the HYP smeared link and the all-to-all propagators. Our purpose is 
two- fold. The primary purpose is to perform the first high precision study of in nf = 2 
unquenched QCD, which serves a reference point for future studies with better control 
over the systematic errors. The secondary goal is to understand in what conditions the 
above methods apply efficiently. We observe the dependence of the statistical errors on the 
time and the numbers of low-lying eigenmodes, as well as their behavior against variation 
of the quark mass and the lattice spacing. This will help us to understand in which region 
of parameters the method can give good control over the statistical errors, which will also 
be useful to precision calculations of other physics parameters for heavy-light systems. 

This paper is organized as follows. In section [21 we describe the method to obtain 
B*Btc coupling from the B meson matrix element. Section E] explains our simulation 
details. In this section we first arrive at our final result for the B*Btt coupling with our 
best parameter setting. Then in Section HI the efficiency of the low mode averaging is 
examined in detail. Conclusion is given in Section [5j 
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2 Lattice observables 



The Lagrangian of heavy meson effective theory is given as 

L = -Tr [Hiv ■ DH] + g b Tr [HHA^ • 7m75 ] + 0(1/M), (7) 

where the low energy constant % is the B*Bir coupling, v is the four- velocity of the 
heavy-light meson B or B*, and H, D^, are described by the B, B* and it fields as 

H = ^(l + lf , Vfl )(iB l5 + B;^), £ = exp(i7r//), (8) 

= d^-^d^ + td^), Ap = (9) 

The B*Btt coupling can be obtained from the form factor at zero recoil which corresponds 
to the matrix element 

(B*(p B *, A)|A|£(p B ))| pV=fe=0 = + m B .)^i(? 2 = 0)ef } , (10) 

where Ai(q 2 = 0) is the matrix element of the transition from B to B* at zero recoil with 
axial current Ai = ip^jiip and A stands for polarization [2]. In the static limit, 

9oo = A^q 2 = 0) (11) 

holds. The matrix element (B*\A I1 \B) at the zero recoil can be obtained from the ratio 
of 3-point and 2-point functions, R(t). 

wm°m = lim mu) , (12) 

2m b t,t A -+ao 

where 

(Oh*(t + t A )A(t A )O B (0)} _ C 3 (t + t A ) 

[,A) (o B (t + t A )o B (o)) -c 2 (t + t A ) [i6) 

with Ob and O b * are some operator having quantum numbers of the B and B* mesons, 
respectively. We apply the smearing technique to enhance the ground state contributions 
to the correlators as 

<D B (t,x) = ^0(f)g(t,£ + f) 7 5 /i(t,f), (14) 

r 

B *(t,x) = ^^r^f + rKM^), (15) 
f 

where 4>(x) is the smearing function. 
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The lattice HQET action in the static limit is defined as 

S = J2 H x )^f^ \ h (x) - U\{x - l)h{x - 4))j , (16) 

x 

where h(x) is the heavy quark field. The static quark propagator is obtained by solving 
the time evolution equation. As is well known, the HQET propagator is very noisy, and it 
becomes increasingly serious as the continuum limit is approached. In order to reduce the 
noise, the Alpha collaboration P, [7] studied the HQET action in which the link variables 
U^{x) are replaced by the smeared links W^(x) in order to suppress the power divergence. 
They found that the noise of the static heavy-light meson is significantly suppressed with 
so-called HYP smearing [8]. 

The statistical error is further suppressed by applying the all-to-all propagator tech- 
nique developed by the TrinLat collaboration [9j. Defining the Hermitian lattice Dirac 
operator Q = 75-D, where D is the lattice Dirac operator, the quark propagator S q (x,y) 
is expressed by the inverse of the Hermitian Dirac operator Q = Q^ 1 as 

S q (x,y) = Q(x,y)~f 5 . (17) 

We divide the light quark propagator into two parts: the low mode part and the high 
mode part. The low mode part can be obtained using low eigenmodes of Hermitian Dirac 
operator Q. The high mode part can be obtained by the standard random noise methods 
with time, color, and spin dilutions. With the projection operators into the low and high 
mode parts, 



JVev 

5% (< >(a;)®v0t(y), P 1 = l-P 



respectively, the propagator can be decomposed into two parts as 

Q = Q + QP U (19) 

Qo(x,y) = J2 T v^(x)(E)v^(y), (20) 
i=i A * 

T r j 

where N r is the number of random noise and j is the index for dilution to label the set of 
time, spin and color sources, j = (t ,ao,a ). The low mode part Qo is constructed from 
the eigenvectors v ^ with their eigenvalues A«, which are to be obtained at a preceding 
stage. As the random noise vector for the high mode part, we adopt the complex Z 2 noise. 
The random noise vector with dilution is given as 

= V[r](x) a J t , to 6 a>ao 6 a ,a - (22) 
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Table 1: The simulation parameters. The values of the lattice spacing and the pion mass 
are from Ref. |13j . 

■0 is given as 

Ar]( x ) = ^2(Q p i)( x iy)vir](y), (23) 

y 

which is obtained by solving a linear equation Qif)[ r ] = PiT][ r ]. Further details of the 
computation methods are given in Ref. [T2] . 

Combining these propagators, we can obtain the 2-point functions for the heavy-light 
meson which are averaged all over the spacetime. Similarly, the 3-point functions can be 
divided into four parts: low-low, low-high, high-low and high-high parts. 

3 Results 

3.1 Simulation setup 

Numerical simulations are carried out on 12 3 x 24 lattices at (3 = 1.80 and 16 3 x 32 
lattices at (3 = 1.95 with two flavors of 0(a)-improved Wilson quarks and Iwasaki gauge 
action. We make use of about 100 to 150 gauge configurations provided by CP-PACS 
collaboration [13] through JLDG (Japan Lattice DataGrid). We use the 0(a)-improved 
Wilson fermion for the light valence quark with the masses set equal to the sea quark 
masses. We use the static quark action with the HYP smeared links with the smearing 
parameter values («i, 02,0:3) = (0.75,0.6,0.3) (HYP1) [HI El- The B and B* meson 
operators are smeared with a function 0(r) = exp (—0.9 f), where f is the distance between 
the heavy quark and the light quark in lattice units. The configurations are fixed to the 
Coulomb gauge. 

We obtain the low-lying eigenmodes of the Hermitian Dirac operator using implicitly 
restarted Lanczos algorithm. The low mode parts of the correlation functions are com- 
puted with N ev = 200 low-lying eigenmodes, except for the case of k = 0.1375 at ft = 1.95 
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Figure 1: The effective mass plot of the 2-point and 3-point functions at /3 — 1.80. Top 
left, top right, bottom left, bottom right panels correspond to k — 0.1409, 0.1430, 0.1445, 
0.1464, respectively. 



which is obtained with N ev = 0. The reason of this choice will be explained in Sec. HI 
The high mode parts of the correlation functions are computed with complex Z2 random 
noise vector with N r = 1. The number of time dilution for each configuration are set to 
N to = 24 at P = 1.80 and iV to = 32 at (5 = 1.95, respectively. This setup is based on the 
experience from the work by Negishi et al. p2] . 

3.2 Correlation function and effective mass 

Figures [1] and [2] show the effective mass plots for the 2-point and 3-point functions . We 
find that the 2-point functions exhibit nice plateaux at t > 4 for (3 = 1.80 and at t > 5 
for (3 = 1.95. From this result we take — 5 for (5 — 1.8 and tA = 6 for (3 — 1.95 as a 
reasonable choice for the time separation between the current Ai and the B meson source. 
We also find that the effective masses of the 3-point functions give consistent values with 
those of the 2-point functions. We fit the 2-point and 3-point functions to exponential 
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Figure 2: The effective mass plot of the 2-point and 3-point functions at /3 — 1.95. Top 
left, top right, bottom left, bottom right panels correspond to k — 0.1375, 0.1390, 0.1400, 
0.1410, respectively. 



functions with single exponent as 

C 2 (t) = Z 2 exp (-E stat t), C 3 (t) = Z 3 exp (-E stat t), (24) 

where Z 2 and Z 3 are constant parameters and E stat corresponds to the energy of the 
heavy-light meson. The fit ranges are chosen appropriately by observing the effective 
mass plots as listed in Table El The bare B*Btt coupling can be obtained by the ratio of 
the fit parameters as g^ re = Z 3 /Z 2 . Alternatively, the B*Bn coupling is also extracted 
from the ratio of the 3-point and 2-point functions, C 3 {t)/C 2 {t), as shown in Figures [31 
We find that the fit of the ratio C 3 (t)/C 2 (t) to a constant value gives consistent value 
with the value of Z 3 /Z 2 . Since the statistical accuracy is better, we employ Z 3 /Z 2 to 
determine in the following analyses. The results are summarized in Table [21 
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Table 2: The numerical results of the heavy-light meson energy aE, the ratio of the 3-point 
and 2-point functions Z3/Z2, and #00. (^min'^max) an d (^min' ^max) are the fit ranges for 
the 2-point and 3-point correlators, respectively. For the values of g^, only the statistical 
errors are quoted. 



3.3 Physical value of the B*Btt coupling and chiral extrapolation 

The physical value of the B*Btt coupling is obtained by multiplying the bare value by 
the renormalization constant. We use the one- loop result of the renormalization factor 
for the axial vector current 

M = 2ku Z a (l + b A ^j Af\ (25) 
1 

( 8412 \ * 
«o=( 1 -^— > &a = 1 + 0.0378^0/), 



ft 



where the gauge coupling g|^(/i) = 3.155, 2.816 and Z A = 0.932, 0.939 for (3 = 1.80 and 
1.95, respectively, as given in Ref. [13]. We arrive at the results of for our k values in 
Table H 

We take the chiral extrapolation of the B*Btt coupling employing the following fit 
functions: 



(a) 


9 {a \ml) 


= g(Q)+A 1 ml, 




(b) 


9 {b \ml) 


= g(0)+A im l + A 2 (ml) 2 , 




(c) 




= 9(0) 


Sir 2 


+ A x m\ + A 2 {mlf 



corresponding respectively to (a) the linear extrapolation, (b) the quadratic extrapolation, 
and (c) the quadratic plus chiral log extrapolation where the log coefficient is determined 
from ChPT |HJ [151 [16] . We use three lightest data points for the fit (a), while all the four 
points for (b) and (c). We obtain physical values of the B*Bn coupling in the chiral limit 
as c/oo = 0.57(1), 0.57(2), 0.52(1) at f3 = 1.80 and goo = 0.548(6), 0.529(10), 0.480(8) at 
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Figure 3: The ratio of the 3-point and 2-point functions versus t at /3 — 1.95. Correspon- 
dence of the panels and the values of k is the same as in Fig. [2j 



(3 = 1.95 from the fit (a), (b), (c), respectively. Figure H] shows these chiral extrapolations. 
We take the average of the results from the linear fit and the quadratic plus chiral log 
fit as our best value and take half the difference as the systematic error from the chiral 
extrapolation: 

9K = 0) = 0.543(5) stat (26) chiral at = 1.80, (26) 
= 0.516(5) sta t(31) chiral 

at {3 = 1.95. (27) 

Since we have only two lattice spacings, naive continuum extrapolation would not give 
a reliable result. However, the results at these two lattice spacings are consistent within 
quoted errorsQ Therefore, we take the result at (3 — 1.95 as our best estimate for the 

1 We do not observe a large scaling violation for as opposed to the case of / w by CP-PACS 
collaboration. A possible explanation is that the large scaling violation for /„. might come from the 
perturbative error of ca, which is the (9(a)-improvement coefficient of the light-light axial vector current. 
If this is the case, f w receives significant systematic error from the 0(a)-improvement term for the 
temporal axial vector current ac^^P due to the chiral enhancement, whereas g m does not receive such 
a systematic error since the corresponding term ac^^P drops out from the matrix element by the sum 
over the space. 
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Figure 4: The chiral extrapolation of the physical B*Btt coupling at (3 — 1.80 (left panel) 
and j3 = 1.95 (right). 



physical value of g^, and estimate the discretization error of 0((aA) 2 ) by order counting 
with A ~ 0.3 GeV. Including the perturbative error of 0(a 2 ) also by order counting, our 
results for is 

#T 2 = 0.516(5) stat (31) chiral (28) pcrt (28) disc . (28) 

In our study, each of the chiral extrapolation error, the perturbative error, and the dis- 
cretization error is about 6% level. The perturbative error can be removed by employing 
the nonperturbative renormalization such as the RI-MOM scheme, which is successfully 
applied to the light-light axial vector current. The discretization error could be reduced 
by computing on finer lattices. For example, an order counting estimate suggests that the 
discretization error would be reduced to about 2% on the configurations of CP-PACS at 
(3 = 2.10. In contrast, it is not straightforward to control the chiral extrapolation error. 
It is definitely necessary to use recent unquenched configurations with smallest pion mass 
m-n ~ 0.3 GeV. More predominant approach is to employ a fermion formulation possess- 
ing the chiral symmetry, such as the overlap fermions, which makes the extrapolation 
theoretically more transparent. 

4 Applicability of the low mode averaging 

In this section, we examine under what condition the all-to-all propagator technique, in 
particular the low mode averaging, is efficient to reduce the statistical error. This would 
give us a guide to extend our computation to unquenched simulations with smaller quark 
masses and finer lattices. We mainly investigate the case of (3 = 1.80 in the following. 
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Figure 5: The low-lying eigenmode distribution for various k at /3 
configurations. 
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Figure 6: The low and high mode contributions to the 2-point correlators versus t at 
(3 = 1.80. The left and right panels show the results at k — 0.1464 and 0.1430, respectively. 



4.1 Observation of the noise to signal ratio 

Figure [5] shows the distribution of about 250 lowest eigenmodes for each k. Since the 
contribution of each mode to the correlator is multiplied by 1/Aj, one naively expects 
that with a fixed number of modes the low mode averaging should be particularly effective 
for the smallest quark mass. This is indeed true for the 2-point and 3-point heavy-light 
meson correlators. Figure [6] represents the low and high mode contributions to the 2-point 
correlators. For the smaller quark mass, the low mode part is indeed more dominant the 
2-point correlator in the earlier t. 

In Figure [7] we show the comparison of the noise to signal ratio of the 2-point functions 
with different number of low eigenmodes, iV ev = 0, 50, 100, 200 for k = 0.1464 and 
k = 0.1430 at /3 — 1.80. For the smallest quark mass, the statistical error of the 2-point 
function is 1.5-2 times improved as N ev is changed from to 200. While this is not a 
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Figure 7: The time dependence of the noise to signal ratio for k = 0.1464 (left panel) 
and k = 0.1430 (right) at (3 = 1.80 with 40 configurations. 



drastic improvement, comparing the costs to determine the low-lying eigenmodes and to 
solve the quark propagator (the latter is 8-10 times larger than the former), there is still 
an advantage to adopt the low-mode averaging. Projecting out the low-lying modes also 
improves the cost of solving the quark propagator. These effects are amplified as going 
to smaller quark mass region. 

In the region of larger light quark mass, however, the situation is different. The right 
panel of Fig. [7] shows the noise to signal ratio of the 2-point function for k = 0.1430. The 
noise to signal ratio for t < 7 achieves about factor 1.3 improvement in the statistical 
error as we change N ev from to 200. For t > 7, on the contrary, the noise to signal ratio 
with N ev starts to grow more rapidly than that with N ev = which keeps to grow 
steadily. As the result, the low- mode averaging deteriorates the statistical accuracy at 
this and larger light quark masses. 

4.2 High mode and low mode contributions to the noise 

To investigate the origin of this behavior, we examine the high mode and low mode 
contributions to the error of correlator separately. 

In general, by projecting out larger number of low modes, the high mode part de- 
creases. Therefore one might naively expect that the error also decreases. This is indeed 
the case for k = 0.1464. However for k = 0.1430 such a naive expectation does not hold. 
Figs. [8] shows the time dependence of the error for the high mode part of the 2 point 
correlator with various values of N ev at k — 0.1464 and k = 0.1430. As is displayed in the 
right panel of Fig. [HI although the error of the high mode contribution to the correlator 
at k — 0.1430 does decrease at t < 7 with larger N ev , the errors of the high mode part at 
t > 7 with N ev = 50, 100, 200 exceed that with N ev = 0. This clearly indicates that for 
t > 7 both the high and low mode parts of the correlator individually have large errors, 
but when they are combined the error of the total correlator becomes small. In this situa- 
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Figure 8: The time dependence of the noise from the high mode part for k = 0.1464 (left 
panel) and k = 0.1430 (right) at (3 — 1.80 with 40 configurations. 
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Figure 9: The noise to signal ratio of the low and high mode parts for the 2-point 
correlators for k = 0.1430 at (3 = 1.80 with 40 configurations. The signal part is always 
taken to be the total correlator using only the noisy estimator. The projection is made 
with 50 eigenmodes. 



tion if the low mode part is improved by the low mode averaging, which reduces the error 
by certain factor, the error of the high mode part of the same size remains unreduced and 
dominates the error of the correlator. 

To see it more clearly, let us decompose the 2-point correlator computed by the noisy 
estimator (corresponding to N ev =0) into the high and low mode parts. The total corre- 
lators are computed only with the noisy estimator (denoted as 'total(noise)'). The high 
and low mode parts ('high(noise)' and 'low(noise)') are separately computed using the ex- 
actly same random source as for the total correlators but projected into the high and low 
mode spaces with the projection operators Pi and P , respectively. Figure [9] displays the 
statistical errors from the low and high mode parts normalized with the total correlator, 
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P= 1.80, K = 0.1464 

• — • Total (noise) 

□ Low mode (N, v =200) 
A — A Low mode (N - 50) 
y— V Low mode (N - 16) 
— ^ Low mode (N - 4) 
X — x Low mode (N - 1 ) 





Figure 10: The noise to signal ratio of the low and high mode parts for the 2-point 
correlators at f3 — 1.80 with 40 configurations. The signal part is always taken to be 
the total correlator using only the noisy estimator. The left and right panels show the 
dependence of the error of 'low(noise)' on N ev for k = 0.1464 and k = 0.1430, respectively. 



in the case of N ev = 50 at k = 0.1430. For comparison we also show the error of the low 
mode part determined with low mode averaging ('low(eigen)'). This figure confirms that 
the fluctuations of the low and high mode parts are almost the same size and compensate 
in the total correlator so as to give much smaller error. The errors of the low mode parts, 
E n oiseif) and E eigen (t), exponentially grow with similar rates, while different from that of 
the total correlator. 

Such a behavior continues as we decrease N ev even down to a few N ev . The right 
panel of Fig. [10] shows the case of N ev = 1, 4, 16, 50, 200 for k = 0.1430, where the 
fluctuations of the low mode part, E noise (t), grows similar rates, while absolute values are 
shifted downward. The left panel of Fig. [TU] shows the case for k = 0.1464. We observe 
that the low mode part ('low(noise)') does not exceed the signal by large amount, which 
explains the behavior observed in Fig. [7J 

The reason why projecting into the low or high mode part provides a drastic en- 
hancement of the error is still unknown. Although the phenomena themselves are quite 
interesting and deserves for further studies, in this paper we restrict ourselves within their 
implication to applicability of the all-to-all propagator technique to the static heavy-light 
system. 

4.3 When is the low mode averaging efficient? 

We have seen that the low mode averaging is efficient only if the error of the noisy 
estimator (not the correlator itself) is dominated by the low mode part. Once the errors 
from the high and low mode parts of the correlator start to exceed the error of the total 
correlator, the low mode averaging is no longer effective but it makes the situation even 
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t[fm] 

Figure 11: The time dependence of the noise to signal ratio of the 2-point correlators at 
(3 = 1.95. The results are determined with N cv = 200, while for k = 0.1375 the N ev = 
result is also displayed. 



much worse. Our result implies that at (3 = 1.80 and k < 0.1430, the rapid growth of the 
error at t > 7 in Fig. [7| signals the breakdown of the above condition. For k < 0.1430, 
the noisy estimator without the low mode averaging works better. Thus the low mode 
averaging is only efficient in the small quark mass region. However, since we have already 
taken the data and they provided satisfactory statistical accuracy of 2% level, we adopted 
the result with the low mode averaging propagator at (3 — 1.80. As for (3 = 1.95, the 
low mode averaging has not provided sufficient statistical accuracy for k = 0.1375. Thus 
we adopted the noisy estimator without the low mode averaging at this k as was already 
noted in the previous section. 

The results at /3 — 1.95 are displayed in Figure [TTJ The figure shows the noise to 
signal ratio of the 2-point correlators at each k against t in physical units. For k = 0.1375 
both the results with N ev = and 200 are shown, and the former indeed exhibits smaller 
statistical error. For all the values of k with N ev = 200, the slopes of the exponential 
growth rate of the noise to signal ratio change around t ~ 16, and beyond that t the slopes 
become steeper as the quark mass increases. This behavior is clearly explained with the 
breakdown mechanism of the low mode averaging mentioned above. 

We can also extract a hint on the lattice spacing dependence of the statistical accuracy 
by comparing f3 = 1.80 and f3 = 1.95. Comparison of Figs. [7] and [11] implies that the 
noise to signal ratio is similar or even smaller for finer lattices. This is partly explained 
by the fact that as going the finer lattices one has the larger number of lattice points (if 
the volume is kept unchanged) which are used for the all-to-all propagator. As observed 
in Figs. [Hand [3 the statistical errors at (3 = 1.80 rapidly increases beyond t ~ 8, which 
corresponds to t ~ 1.6 fm in physical units. Thus the low mode averaging breaks down 
almost at the same physical distances at these two lattice spacings. This implies that also 
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Figure 12: Comparison of g^ with other calculations [3], H], [12] . In our results, the small 
and large errors represent the statistical error and the perturbative error, respectively. 



on finer lattices of a ~ 0.1 fm, one can extract precise values of B*Bn coupling from the 
region t < 1.6 fm by applying the same methods as this work, while careful tuning of the 
smearing function would be indispensable. 



5 Conclusion 

In this paper, we computed the B*Bir coupling on unquenched lattices using the HYP 
smearing and the all-to-all propagators. Using the low mode averaging with 200 eigen- 
modes, the statistical errors are kept sufficiently small for smaller quark masses. On the 
other hand, as was investigated in Sec. H]in detail, the low mode averaging is not efficient 
for larger light quark mass region, where the simple noisy estimator provides better pre- 
cision. In either case, the statistical error is controlled below 2% level in the chiral limit. 
We obtained consistent results at two lattice spacings. Our best estimate of the B*Bn 
coupling in the static limit is represented in Eq. ( 1281) . Figure fT2l compares our results 
with other recent works on the B*Bit coupling [2JIHEI2]. The improvement in statistical 
precision is drastic, which proves the power of the improvement techniques employed in 
this paper. 

For future prospects, better control over the systematic error from the chiral extrap- 
olation is indispensable. For this purpose, the configurations with dynamical overlap 
fermions by JLQCD collaboration would be a good choice [HI [EHl [20] . In order to obtain 
cjb at the physical bottom quark mass, one needs to understand the mass dependence of g. 
Simulations with the charm quark mass region and interpolation with the static limit are 
desired. The methods adopted in this work are in principle also applicable to other weak 
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matrix elements of the B mesons, such as Bb, and the form factors, and expected to 
provide high precision results required in precision flavor physics. 
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